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We observe multiple stable states of nuclear polarization and nuclear self-tuning over a large range 
of fields in a double quantum dot under conditions of electron spin resonance. The observations can 
be understood within an elaborated theoretical rate equation model for the polarization in each of 
the dots, in the limit of strong driving. This model also captures unusual features of the data, such 
as fast switching and a 'wrong' sign of polarization. The results reported enable applications of this 
polarization effect, including accurate manipulation and control of nuclear fields. 



Great experimental progress in the last decade en- 
abled the confinement, initialization and read-out of sin- 
gle spins in quantum dots [1 . Controlled coherent single- 
spin rotations — a key ingredient for quantum manipu- 
lation — were demonstrated recently using the electron 
spin resonance (ESR) [31 |U IS] . The weak hyperfine 
coupling of the electron spin to the nuclear spins in the 
host material appeared to be of great importance in this 
field. It was identified as the main source of qubit de- 
coherence and provides a significant hybridization of the 
spin states [TJIH]. This has stimulated intensive theoret- 
ical and experimental research focusing on nuclear spin 
dynamics in quantum dots [51 nilllTTl [HI [131 [Ti] . 

Overhauser pointed out already in the the 1950s [T5j 
that ESR may provide the buildup of significant nuclear 
spin polarization. Indeed, most ESR experiments on 
quantum dots, aimed at demonstrating electron spin ro- 
tations, also clearly demonstrated dynamical nuclear spin 
polarization (DNSP) [21 [31 SI [5] . 

For ESR driving of a single spin in an almost isolated 
quantum dot, or an ensemble of such dots, the scenario is 
similar to that of the usual Overhauser effect: the direc- 
tion of DNSP is parallel to the spin of the excited elec- 
trons [15l[16]. Recent ESR experiments on self-assembled 
quantum dots have confirmed this picture [5] , and a sim- 
ilar reasoning holds for spin experiments with optically 
pumped dots |17j . In some cases, a bistability has been 
observed: Under the same conditions, the nuclear spins 
in the dot can be either polarized or unpolarized [111 . 

However, several issues can complicate the situa- 
tion. In recent ESR experiments in double quantum 
dots [11 [3l [4] (i) electrons participate in transport dur- 
ing ESR driving, and (ii) there can be different nuclear 
spin dynamics in the two dots. Furthermore, a driving 
magnetic field is in practice accompanied by an electric 
field which modulates the electron-nuclear spin coupling 
at the resonance frequency [Ij . All this makes a straight- 
forward extension of existing models \1& impossible and 
promises richer and more interesting physics, which we 
indeed reveal. 

In this Letter, we report a study of ESR in a dou- 
ble quantum dot focusing on DNSP. We have observed 
multiple stable states of nuclear polarization (up to four 
states), not seen in single-dot experiments, nuclear self- 



tuning to the ESR condition over a large range of mag- 
netic fields (> 100 mT), and a sign of DNSP opposite to 
that following from the Overhauser argument. We iden- 
tify the most probable mechanism governing DNSP and 
present a theoretical model explaining our findings. The 
results reported enable applications of this self-tuning ef- 
fect, including accurate manipulation and control of the 
nuclear polarization [18J and use of this for improving 
the electron spin coherence time, possibly by orders of 
magnitude. 

The double quantum dot system is electrostatically 
defined in a two-dimensional electron gas, located 90 
nm below the surface of a GaAs/AlGaAs heterostruc- 
ture, by applying negative voltages to metal surface 
gates. The dots are tuned to the Pauli spin blockade 
regime [13] , where the transport sequence of charge states 
is (1,1) (0,2) (0,1) (1,1), (n,m) denoting the 
charge state with n(m) excess electrons in the left(right) 
dot. The current through the double dot depends on the 
spin orientation of the electrons in the (1,1) state since 
the only accessible (0, 2) state is a spin sing let (Fig.[l]). 

Magnetic spin resonance is achieved by sending an 
alternating current through a coplanar stripline (CPS) 




FIG. 1: Double dot setup, (a) The double quantum dot is cou- 
pled to two leads. Due to a voltage bias, electrons can only 
run from the left to the right lead, implementing the transport 
sequence (1,1) (0,2) (0,1) (1,1). (b) Energy dia- 
gram. The four possible (1, 1) states differ in spin projections 
on the quantization axes (red arrows). Under ESR conditions 
the axes can be different in the two dots and do not coincide 
with the direction of the external magnetic field. These states 
are coherently coupled (green arrow) to the (0, 2) singlet that 
decays quickly (broadened line), leaving the system in (0, 1). 
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FIG. 2: (a) Magnetic field sweeps for uj fixed at 350 MIIz. Upper panel: Magnetic field sweep from low to high values resulting in 
an ESR peak width exceeding 100 mT. Lower panel: Sweep in the opposite direction, showing a much narrower ESR peak |20) . 
The nominal resonance condition B^es = fi-i^/gfJ'B is met at Bq m 71 mT for u = 350 MHz and g = 0.35 [2] (see dashed line). 
Note that in both traces the nuclear bath is unpolarized at the onset of electron spin resonance [21]. (b) Multiple values of 
the current through the double dot approximately at resonance. The current switches between at least three stable values on 
a time scale of seconds to minutes. The three panels correspond to three different values of the energy level detuning Alr 
(increasing from the bottom to the upper panel). The values given for Alr may have a constant offset, as photon assisted 
tunneling processes broaden the interdot transition which makes it difficult to separate resonant and inelastic transport. In 
both (a) and (b) the lowest value of current was subtracted as offset. The data in (b) were taken for a larger Fin and Fout than 
the data in (a). 



which lies on top of the surface gates, separated by a 
thin dielectric layer. This current produces a small os- 
cillating magnetic field _Bi ~ 1 mT perpendicular to the 
external magnetic field -Bq — 100 niT. The experimental 
data are obtained with the same device and in the same 
measurement run as the data presented in Ref . [2 . The 
difference is that the device is tuned to a higher interdot 
tunnel coupling and coupling to the right lead. 

When we apply a continuous wave RF current with 
fixed frequency uj to the CPS and sweep the exter- 
nal magnetic field Bq passing the resonance condition 
Bq = Bj-cs = huj/gfiB, we make a remarkable observa- 
tion. One would expect that the resonance manifests 
itself as a peak in the current |2]. Indeed, if the external 
field is swept from low to high values, the current jumps 
up upon achieving the resonance condition. Unexpect- 
edly, this resonant response extends over a wide range of 
magnetic fields, that exceeds i?rcs by a factor of 2 (see 
Fig. |2|l upper panel). If the field is swept in opposite 
direction (Fig. [2^ lower panel), the current remains low 
till Bq is several mT above B^^s-- This indicates a strong 
hysteresis for Bq > i?rcs, whereas the hysteresis below 
-Bics is much less pronounced. 

Another unexpected observation is made at fixed Bq « 
Bros- Instead of a single value of the current correspond- 
ing to the maximum value of the ESR satellite peak, we 
observe clearly distinguishable multiple stable values of 
the current. Switching between these values gives rise to 
a random telegraph signal (RTS) with time scales rang- 
ing from seconds to minutes. Typical time-resolved mea- 
surements of the RTS are presented in Fig. for three 
different values of the energy level detuning Alr (Fig.jlJ. 

We associate both the hysteresis and RTS with DNSP 
induced by the non-equilibrium electron spin dynamics 



under conditions of ESR and transport in the dots. Nu- 
clear polarization is known to provide an extra effec- 
tive magnetic field Bpf acting on the electron spin |15j . 
Where high current is observed in the hysteresis re- 
gion, this extra field should be such that the total field 
Bq + Bn ~ -Bros, i-e. the nuclear field 'tunes' the system 
to the resonance condition jl6j . Low current indicates 
that the total field Bq + B]\[ significantly deviates from 
Bj-cs'- The nuclei are unpolarized. Both polarized and 
unpolarized states are stable in the interval of hystere- 
sis. Fluctuations of any kind could provide spontaneous 
switching between stable states, leading to the RTS. 

A number of experimental details does not fit into this 
simple picture. Firstly, there are multiple values of the 
current observed, three are clearly visible in Fig. [2|d (la- 
beled A-C). This implies multiple stable states of nuclear 
polarization with a total field close to i?ros- Actually, we 
think that the RTS traces provide evidence for the ex- 
istence of a fourth state. There is a number of current 
dips observed (labeled D) too big to be statistical fluctu- 
ations. We interpret those dips as signatures of a fourth 
state that decays on the scale of a second, i.e. different 
from state A, which decays on a larger time scale. Sec- 
ondly, switching between the different current levels is 
rather fast. The nuclear spin dynamics are known to be 
slow, with a typical relaxation time t„ ~ 15 s [7l[T2|[T3]. 
If the current is a direct measure of the nuclear polariza- 
tion, then why is the duration of the switching events so 
short? 

A third point is the sign of the polarization. Usually, in 
ESR experiments the dominating mechanism of DNSP is 
described by the Overhauser effect: The ESR excitation 
drives the electron spin(s) out of equilibrium, and hyper- 
fine induced electron-nuclear spin exchange is one of the 
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mechanisms contributing to electron spin relaxation. As 
reasoned by Overhauser, on grounds of spin conservation, 
the direction of nuclear polarization should be parallel to 
the spin of excited electrons, whatever its orientation is 
with respect to the magnetic field applied. This is the 
case for most DNSP experiments, e.g. [3J |5J [TT]. Given 
the negative (/-factor and positive hyperfine coupling in 
GaAs |22], this would give a Bm parallel to Bq [16 . In 
our experiment, its direction is clearly opposite, as high 
current is seen for Bq > JSres- AH three points are cap- 
tured by the theory given below. 

The electron spin S and nuclear spins in each dot 
are coupled by hyperfine interaction |22j 

{2^'4' + s+ik + s-i^} , (1) 

k 

where the sum runs over all TV 10^ nuclei in the dot. 
The energy Ak is proportional to the probability to find 
the electron at the position of nucleus k, Ak — 10~^° eV. 
With an external field applied in the z-direction, the 
'fiip-fiop' terms S^I^^ provide spin exchange between 
the electrons and nuclei. Owing to energy conserva- 
tion, these exchange transitions must be second-order 
processes involving a mechanism supplying or absorbing 
the excess Zccman energy. Conventionally, the electron- 
nuclear spin exchange is due to the time-independent 
hyperfine coupling A^- However, as recently has been 
pointed out [2, J;^, in this setup a significant a.c. elec- 
tric field moves the electrons in the dots with respect 
to the nuclei. This can be accounted for by introducing 
a time- dependent component in the hyperfine coupling 
Ak Ak + ifce'"* + A^.e^*"*. We estimate that under 
the present conditions Ak/Ak — 0.1 |4l 123]. 

We have considered six candidate mechanisms for 
DNSP [21], assuming a saturated ESR. We concluded 
that the dominant one involves the time-dependent hy- 
perfine coupling, which allows for 'photon assisted flip- 
flops'. These flip-flops not have a preferred direction set 
by a large energy mismatch: now the spin asymmetry is 
now provided by internal spin relaxation causing the spin 
ground state (parallel to the external field) to be more 
populated than the excited state. 

The theoretical consideration includes the following 
steps: (i) We consider the four (1, 1) states using a rotat- 
ing wave approximation, assuming a saturated ESR and 
a negligible exchange splitting, i.e. min{t, i^/ALR,} ^ 
Bi,Bn. The eigenstates in a rotating frame are mix- 
tures of spin-up and spin-down states, with a mixing 
angle O^^if — ^ arctan{_Bi_fl/2/i_fl} which can be dif- 
ferent in both dots (see Fig. [TJ, due to e.g. different cou- 
pling of the electrons to the CPS. The Rabi frequency 
in each dot B^ji = giiBB\ ' /h gives the width of the 
saturated resonance, and the ESR frequency mismatch 
fL,R = \91^b{Bq + B^'^)/'h\ — Lo depends on the nu- 
clear polarization in each dot. (ii) We evaluate the tran- 



sition rates between these states to obtain their quasi- 
stationary population and the current through the dou- 
ble dot. We include tunneling (characterized by = 
t^/Fout — 1 — 10 MHz) and single electron spin relax- 
ation |16j (cx Fj. ~ 1 MHz at zero temperature, which will 
be enhanced by a thermal factor UbT / g^sBo = ^ ~ 5, 
in accordance with a lower bound estimate set by the 
typical leakage current of 100 fA). This approach is valid 
in the limit B ^ F^ (iii) We compute the rates of 
hyperfine-induced spin exchange. In the first approxi- 
mation we find rates symmetric with respect to nuclear 
spin, their scale set by Fa ~ A^/(64?i^^Fr) - 0.5 Hz. Be- 
ing symmetric, these rates do not contribute to DNSP. 
They merely enhance the relaxation of the nuclear fields, 
(iv) The small spin-asymmetric part of these rates Fi ~ 
\{Ak/?>TiBf{Vs/0 ^ 10"^ Hz, due to electron spin re- 
laxation, introduces a preferential direction of nuclear 
spin pumping in each dot. (v) We construct equations 
of motion for the effective nuclear fields B^'^ and an- 
alyze the stable states of nuclear polarization given by 
dB^'"/dt = 0. (vi) We use a Fokker-Planck equation to 
give a qualitative analysis of fiuctuations of nuclear po- 
larization and switching rates between the stable states. 
The evolution equation for B^ thus found reads 

^ = -T,B,,P{dL,B.) - { ^ + ^2RiOLM)^B^, (2) 

and the equation for B^ is obtained by permutation of 
L and R. The field Bqv is the Overhauser field of full 
polarization, _Bov ~ 5 T for GaAs. The functions P and 
R are dimensionless functions giving the functional de- 
pendence of the resonant nuclear spin pumping (P) and 
resonantly enhanced nuclear spin relaxation (i?) on the 




FIG. 3: (a) Time-derivative dB^/dt at the edge of the hys- 
teresis interval -Bo ~ Bres (green) and in the middle of the 
interval Bo « -Bres -I- 0.5 |B™| (blue). (b,c) Close-up at 
resonance. The curves consist of the usual relaxation (lin- 
ear slope) which is resonantly enhanced (dashed lines), and 
spin pumping that adds a two-peak shape near the resonance. 
The circles indicate the stable states of nuclear polarization. 
We used Fi/Ti = 0.043, FaT^ ^ 5, 6r ^ 0, ^Tr/Ts = 0.75, 
and assumed equally strong coupling Ak of all nuclei to the 
electron. 
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FIG. 4: Stable polarizations in the plane (B^, B§), for the cases (a) Bo — Bros ~ Bi and (d) Bo — Bros ~ 0.5 jB]^"'"!. A contour 
plot of the current is included, the gray shade indicating the region with highest current. Switching between the stable points 
gives rise to RTS as presented in (b) and (c). A qualitative difference is that the point e in (d) is 'isolated', i.e. having switched 
to e, the system will never switch back. In (a) an asymmetry in Bl,_b and N^^ji is included, resulting in four different current 
levels for a-d, whereas (d) is plotted assuming a symmetric double dot. Note the different scales at the axes in (a) and (d). 
The same plots (a) and (d) can be found in the Supplementary Material where we included the local nuclear spin dynamics as 
a vector field. 



mixing angles and on Fs/^Fr, and have a maximum ~ 1. 
While R is roughly Lorentzian- shaped, the function P is 
zero far from resonance 9 — > {0,7r}, reaches maximum 
at I /I ~ B, and falls off to zero again at the resonance 
9 = 7r/2. This resonant dip is due to the vanishing of 
electron spin polarization at the saturated resonance. In 
Eq. ([2]), the terms proportional to —B^ give nuclear spin 
relaxation: The first term presents the usual r„ while the 
second term gives a resonant enhancement owing to spin 
exchange with electrons. Nuclear spin pumping is given 
by TiBovP (~ 50 mT/s, much faster than the sweep rate 
in Fig.|2^), with a sign opposite to that following from the 
Overhauser reasoning: Spin exchange under conditions 
of electron transport is mostly due to electrons polarized 
along the direction of the external field. The shape of a 
typical pumping curve is shown in Fig. |3] 

We are now also able to understand the extended in- 
terval of hysteresis: ESR response can be observed as 
long as there exist stable solutions of dBjq/dt = close 
to resonance. Eq. ^ determines the interval of hystere- 
sis as Bros ^ ^0 < -Sres + I^Tv^'^l, where the maximal 
nuclear field is B]^^'' = -Bovri/(F2 + r^^). Using the 
parameters as estimated above we find that F2T„ ^ 10. 

It is the two-peak shape of the pumping curve that is 
responsible for the multiple stable states of nuclear po- 
larization, even at the edge of the hysteresis interval. If 
Bq « Bres (Fig. [3] green curve), there are four stable 
states for the double dot system. This is represented 
in Fig. |4^, where the circles indicate the stable points 
in the plane {Bj^,B^). It is now clear how, even close 
to Bq = Bj-cs, the system can have four stable states 
with different current. A rough estimate for the dura- 
tion of the switching between those states is the typical 
distance over the local speed of the spin dynam- 

ics (~ FiBov), giving ~ 10~^ s, which explains the fast 
switching. A typical time trace in this case will look like 



Fig. |4Jd, which is to be compared with Fig. [2]d. 

When increasing Bq, both dots will develop a separate 
third unpolarized stable state (Fig. [s] blue curve), giv- 
ing as many as nine stable points, as presented in Fig. 
|4jj. At higher fields the unpolarized state (labeled e) will 
become isolated from the other stable states: If the sys- 
tem switches to e, it will never switch back (see Fig. |4|:). 
This also has been observed in experiment [H]. When 
subsequently sweeping back from high to low field, the 
barrier for switching back from e to a high-current state 
is again gradually lowered. When the typical switching 
time becomes comparable to the time scale of the sweep, 
one can expect the current to switch to a high value (Fig. 
|2^, lower panel). 

From Eq. ^ we construct a two dimensional Fokker- 
Planck equation to study the stochastic properties of 
the polarizations in more detail |16j. Importantly, due 
to the accelerated dynamics, the fluctuations around 
all polarized states are suppressed as {{ABn)'^) /il'^ w 
(Bi/|B]^^'^|), = {Ak/g^lBfN being the field variance 
in the unpolarized state. Using Kramers' method [24^ we 
derive an expression for the switching rates between the 
stable states. All rates have the exponential dependence 
Fsw oc eyiY>{—aBi\B^^^\l^^}, where a is a numerical fac- 
tor: The rates are suppressed exponentially with a power 
^ B\/{{ABnY) ^ 1. This exponential dependence ex- 
plains the large RTS time scale as well as the strong 
variation with A^j^ in Fig. |2] We calculated the expo- 
nent explicitly for Fgw from a to d in Fig. |4]i. We used 
V2Tn = 10, Vs/£Tr = , and Bo-Bres = 0.5 \B'^^^\ and 
found that a k, 0.72. 

To conclude, we have observed multiple nuclear polar- 
ization states and locking of the ESR condition over a 
large range of magnetic fields in a double quantum dot 
under ESR. We presented a theoretical model that cap- 
tures the existence of these phenomena and their unusual 
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features as fast switching and a 'wrong' sign of DNSP. We 
acknowledge useful discussions with M. Laforest. This 
work was supported by the Dutch Foundation for Fun- 
damental Research on Matter (FOM). 
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SAMPLE 

The experimental data presented are obtained with the same sample as used in reference [T]. A device with the 
same gate pattern as used in the experiment is shown in Fig. |5|l. The two coupled semiconductor quantum dots are 
defined by surface gates (Fig.jsji) on top of a two-dimensional electron gas (2DEG). The oscillating magnetic field that 
drives the spin transitions is generated by applying a radio-frequency (RF) signal generated by a Rohde & Schwarz 
SMR40 source to an on-chip coplanar stripline (CPS) which is terminated in a narrow wire, positioned near the dots 
and separated from the surface gates by a 100-nm-thick dielectric (Fig. [sja). The current through the wire generates 
an oscillating magnetic field Bi at the dots, perpendicular to the static external field i?o and slightly stronger in the 
left dot than in the right dot. 
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FIG. 5: ESR device, (a) Scanning electron microscope (SEM) image of a device with the same gate pattern as used in the 
experiment. The Ti/Au gates are deposited on top of a GaAs/AIGaAs heterostructure containing a two-dimensional electron 
gas 90 nm below the surface. White arrows indicate current flow through the two coupled dots (dotted circles). The directions 
of the external magnetic field and the ac magnetic field are indicated, (b) SEM image of a device similar to the one used in 
the experiment. The termination of the coplanar stripline is visible on top of the gates. The gold stripline has a thickness of 
400 nm and is designed to have a 50 SI characteristic impedance, Zo, up to the shorted termination. It is separated from the 
gate electrodes by a 100-nm-thick dielectric (Calixerene). 

The GaAs/AlGaAs heterostructure from which the samples were made was purchased from Sumitomo Electric. 
The 2DEG has a mobility of 185 x 10^ cm^/Vs at 77 K, and an electron density of 4-5 x 10^^ cm~^, measured at 
30 mK with a different device than used in the experiment. 

Background charge fiuctuations made the quantum dot behavior excessively irregular. The charge stability of the 
dot was improved considerably in two ways. First, the gates were biased by -1-0.5 V relative to the 2DEG during the 
device cool-down. Next, after the device had reached base temperature, the reference of the voltage sources and I/V 
converter (connected to the gates and the 2DEG) were biased by +2 V. This is equivalent to a —2 V bias of both 
branches of the CPS, which therefore (like a gate) reduces the 2DEG density under the CPS. 

The measurements were performed in a Oxford Instruments Kelvinox 400 HA dilution refrigerator operating at a 
base temperature of 35-40 mK. 

MEASUREMENTS 

In both traces in Fig. 2a in the main text, the nuclear bath is unpolarized at the onset of electron spin resonance. 
In the case of the upper panel, Bq is swept just before the measurement from 300 mT to —20 mT in about 40 s. 
During this sweep the nuclear field relaxes (typical relaxation time t„ ~ 10 s) or is even actively depolarized. Residual 
polarization would be indicated by a shift of the zero field peak and the onset of ESR response to a nominal non- 
resonant magnetic field, which is both not observed. In the case of the lower panel the magnetic field is swept from 
low to high magnetic field (to 300 mT) just before recording the trace. In this case there could be polarization still 
present at the beginning of the trace, however in that case that polarization relaxes much faster than the sweep rate 
of 60 mT/min, such that when reaching the resonant field the nuclear spin bath is equilibrated. 
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CANDIDATE MECHANISMS 



Here we describe how we identify the dominating process of hyperfine induced nuclear spin flips. The 'flip-flop' 
terms S^I^ in the hyperfine Hamiltonian are responsible for the exchange of spin between the electron and the nuclei. 
However, as the nuclear Zeeman splitting is 3 or 4 orders smaller than the electron splitting [2\, the states coupled 
by S^I^ are roughly guBBo apart in energy. Therefore, spin exchange is only allowed in a second-order process in 
which some other mechanism supplies or absorbs the excess Zeeman energy. 

This energy difference may (i) be dissipated by an environment f3], or (ii) be given to an electron tunneling out of 
the dot. In case (i) the environment, at sufficiently low temperatures, can only absorb energy, so that the electron 
Zeeman energy can only be reduced. This results in the same sign of DNSP as with the usual Overhauser effect [1]. 
In case (ii), owing to a voltage bias much larger than the Zeeman energy, the change of energy in the course of a 
spin-flip can be of either sign. In this case, a preferential direction of DNSP will be determined by some other spin 
asymmetry of the system. Such an asymmetry may arise from either (ii.a) a difference in spin-flip rates for different 
spin directions (e.g. due to different overlap between initial and flnal states), or (ii.b) different populations of the 
states with different spin directions (e.g. due to internal relaxation processes or differing decay rates). 

Apart from these three mechanisms, there are two more options to choose between: As mentioned in the main text, 
the a.c. electrical component of the exciting field Bi moves the electrons in the dots with respect to the nuclei, and 
this we account for by introducing a time-dependent component in the hyperfine coupling. The time-dependent and 
time-independent couplings will give rise to different flip rates, so this gives us in total six candidate mechanisms. 

Let us first decide on the relative contributions of the time-dependent and time-independent hyperfine couplings, 
Ak and Ak- We compare the strength of second order transition rates, in both cases proportional to the couphng 
amplitude square and inversely proportional to the energy square of the virtual state. While for the time-independent 
coupling this energy is the electron Zeeman energy Ez = gfisBo, it is a much smaller energy for the resonant time- 
dependent coupling, involving the a.c. resonant magnetic field Ej-f = q^bBi. Therefore we have to compare the 
factors {Ak/EzY and {Ak/E^fY. We estimate that for our conditions A^/Ak ~ 0.1 and E^f/Ez = Bi/Bq ~ 0.01, 
and conclude that the time-dependent coupling dominates. 
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FIG. 6: Overview of all mechanisms considered and the corresponding estimates for the nuclear spin pumping rate. The two 
key ingredients for spin pumping are dissipation of the energy mismatch and an asymmetry in spin giving a preferred direction 
for nuclear spin flips. Furthermore, we considered both the effect of Ak and Ak, i.e. time-dependent and time-independent 
hyperfine coupling. We conclude that, under the present experimental conditions, the dominant mechanism involves time- 
dependent hyperfine coupling, energy dissipation by electron transport and internal spin relaxation causing an asymmetry in 
the populations of the electron spin states. 
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To decide upon the other three options mentioned above, we have to compare the spin exchange rates involving 
electron tunneling, characterized by the broadening of the (0, 2) singlet TiFout and the typical decay rate of the (1, 1) 
singlet Ts — f^/Tout, and internal spin relaxation within the dots, characterized by a rate Tr[s], e being the energy 
dissipated. For mechanism (i) we find the scale {Ai^/ E,.f)'^Tr[E,-f], i.e. spin relaxation dissipates the remaining energy 
difference ~ i?rf. In case (ii) the energy is dissipated during tunneling, which takes place with a rate ~ F^, giving 
a scale for the nuclear spin flip rate of {Ak/Ei-i)'^T s- This rate however is symmetric in spin direction, so to find a 
preferred direction of DNSP we need to include an asymmetry: (ii.a) The states are split by ~ Ej-f, so decay to the 
broadened (0,2) singlet introduces a relative difference of ~ E,-{/hTout in the rates, setting the scale of the DNSP 
rate ~ A^F^/i^rfTiFout- (ii-b) Internal spin relaxation competes with tunneling processes, causing an asymmetry in 
the population probabilities of the states of ~ F^ [Ez] /Fg resulting in {Ak/E,f)^Tr[Ez] for DNSP. In Fig. g we give 
a schematic representation of these considerations. We show all mechanisms investigated and give the corresponding 
estimates of the scale of nuclear spin pumping. 

We estimate Ak - 10"" eV, E,f ~ IQ-^ eV, F^ ~ 10 MHz, ?lFout lO""* eV and Tr[Ez] - 100 x Tr[E,{] - 1 MHz, 
resulting in the estimates for the scales of DNSP rate (i) 10~^ Hz, (ii.a) 10~* Hz, and (ii.b) 1 Hz. Based on 
this argument we conclude that mechanism (ii.b) dominates: Electric field assisted hyperfine flip-flops involve the 
absorption and emission of photons with energy tiuj. Close to resonance this effectively reduces the energy mismatch 
of the states involved in a flip-flop from gfisBg to the energy scale of the ESR driving g^sBi- Since this energy 
mismatch is too small to result in a significant nuclear spin pumping rate based on a standard Overhauser argument, 
another spin asymmetry is needed. Internal electron spin relaxation provides this asymmetry: it causes the electron 
spin ground state to be (slightly) more populated than the excited state. This difference in populations combined 
with photon assisted hyperfine flip-flops (which do not have a preferred direction) results in DNSP parallel to the spin 
of the electron ground state. 

THEORY 

Here we will elaborate further on the six steps of the theoretical consideration as sketched in the main text. 

(i) The Hamiltonian for the electron spin operators S^jj in the rotating wave approximation reads 

H = -hfJl - hfnS],, + J [BlSI + BnSf,) , (3) 

L{R) referring to the left(right) dot. The rotating wave approximation is justified by B, \ f\ <C uj- The eigenstates of 
H form the basis {|-f-)^ , \-) l} {\+) ji , with |-|-) = cos 6* |t) + sin 6* |i) and |-) = sin 6* |t) - cos 6* where the 

mixing angle is 9l,r. = 5 a,rcta,n{BL,R/2fL,R}. 

(ii) The master equation includes the decay and relaxation rates, and is justified if B by far exceeds these rates [5]. 

The rates depend on the wave functions of the states involved. Any basis state \n) € , |-l — ) , | — h) , | )} 

decays via the (0, 2) singlet to (0, 1) with a rate F" = |('S'|'^)| Fs, with \S) being (1, 1) singlet. Such a decay process is 
followed by a charge transfer in the left junction (0, 1) |m), whereby all four basis states |m) are re-initialized with 
equal rates Fs/4. Internal relaxation processes are due to coupling to an environment and involve energy dissipation of 
±Ez [3]. We believe that the environment are mainly the electrons in the leads. Their temperature is typically large, 
^ = ksT/Ez ~ 5 1, so we need to consider both emission and absorption rates. They read Fabs — n'B{Ez)Tr[Ez] 
and Fein ~ Fabs + F,.[iJ2], with nsis) being the Bose distribution and Fr[-E2] being the emission rate at zero 
temperature. In the high-temperature limit we find the transition rates FJ?^™ ~ {"^ ^ J2l r I ("^1 ^lr I") P}ri-[^2]- 
We are now able to construct a master equation 

= -F>„ ^TP.n + {F™^>™ - Fr ™P« } , (4) 

and solve it for the quasi- stationary populations pn- These populations gain, via the rates and FJ?"^"^, a resonant 
dependence on /lh on the scale f :^ B and therefore also depend on the nuclear polarizations Bj^' . From the 
populations p„ we can calculate the current through the double dot as /dot = eF^ \ {S\n)\^Pn- 

(iii-v) The rates of electron-nuclear spin exchange are calculated using second order perturbation theory. The 
positive and negative spin flip rates per nucleus in the left(right) dot read [B] 



(5|m) (m|S'^(^^ \n) 



En — En 



Pn- 
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FIG. 7: The current through the double quantum dot /dot as a function of /l and /r. We find high current when only one of 
the two dots is on resonance and low current in the rest of the plane. To generate this plot we used ^ = 5 and Ts = 20rr[-Ez] 
and we subtracted the leakage current far away from both resonances. 



Non-zero diagonal matrix elements such as (+-1-1 will give rise to very small denominators in of the 

order of the nuclear Zeeman energy. Therefore, we have to investigate the contribution of these, possibly dominating, 
terms in another way. We write the second order perturbation in the hyperfine Hamiltonian, 



dp 



H^{t) + H'j,{t), H^{t') + H'j,{t'),p dt', (6) 



where the perturbation is H^j^^{t) ~ i Sfc ^fe{'^L(fl)(^)^irL(K)(*) + ^L{R)(^)^k l{r)(^)}- After separating the time 
scales of the electronic and nuclear spin dynamics, assuming that we can separate the electronic and nuclear part of 
the density matrix as p = pei 'Si Pnuc, and tracing over the electron part of the density matrix, wc find that we can 
write for the time-evolution of the nuclear field in one of the dots 



dBN _ Al (5 



(X"^ - xnBo. - (i?"" + R^'')Bn , (7) 



dt sr? 

where x"'' = -i j\S''{t)S^{t') - S''{t')S''{t))dt' and R"'' = j\S''{t)S''{t') + S''{t')S''{t))dt' , i.e. the susceptibility 
and zero-frequency fluctuations of the electron spin in the dot under consideration. In Eq. ([t]) we left out the 
contributions proportional to the polarization in the x- and y-direction while they are averaged out to zero. We focus 
on the contributions of the diagonal matrix elements and find that ~ X^^ = and i?^^ -I- is only non-zero 
close to resonance, resulting in a resonant enhancement of nuclear spin relaxation. 
We can combine Eqs ([s]) and ^ in an evolution equation for the nuclear fields 

^ = ~T,B,^P{eL,R) - {f + ^2RiOL.R)^ B^, (8) 

where we added a term describing diffusive spin relaxation cx 1 /r„ . The equation for B^ is obtained by permutation 
of L and R. The scales Fi 2 are the same as defined in the main text, i.e. Fi = '^{Ak/8'hl3)'^{Ts/^) and F2 ~ 
A^/(64?i^^F,.), and the functions P and R read 



_ 32a(H- 4a) cos 6Il cos 6'i{COs(6lL - 6iij)sin^ + 4q:{1 + 8a + cos^(6'l - 6';?)}sin^ 26*1, 



^ 16a sin' 0L cos^O l ^ r) + I + Sa{'3 + I6a) 
^ ^'-"^ l + 16a ■ sm^{0L-eii) + l6ail+Aa) ' ^ ' 

where a is the dimensionless variable a = STr[Ez]/^s- For this representation of P we used the high temperature 
limit, i.e. ^ ^ 1 and assumed for simplicity all electron-nuclear spin couplings Ak equal. 

(vi) We also investigated both the switching rates between the different stable states, and the small fluctuations 
near these states. To estimate the fluctuations, we use a two-dimensional Fokker-Planck equation for the distribution 
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FIG. 8: On the edge of the hysteresis interval, Bo — Bres ~ Bi, assuming an asymmetric double dot (Fig. 4a in the main text). 
Here we add vector plots of {B^,B§} in the plane {Bj^,B^) and the current /dot as color background. The circles indicate 
the stable points of nuclear polarization, (a) Overview of the whole region where stable points are expected: Four stable points 
can be distinguished. (b,c) Close ups around the two stable points with high current (corresponding respectively to points c 
and b in Fig. 4a in the main text). From the background colors we can see that the difference in current is ~ 15 %. To generate 
these plots we used for both dots £.Tr-[Ez]/T^ = 0.25. In the left dot Fi/Fz = 3.6 • 10"^, Far^ = 9.04, and Bq - lu = 3.3 Bf, 
and in the right dot Fi/Fa = 16 • 10"^ Far^ = 20.3, and Bo - uj = 7.1 Sf . 



function of the nuclear fields V{Bj^, B§), where —1 < B^'^ /Bqv < 1. To derive the equation, we regard the nuclear 
dynamics in both dots as a random walk on a discrete set of spin values n = — N^), where A^kd is the number 

of nuclei with spin up(down). The DNSP rate Fi only causes transitions from n to n + 1, while the spin relaxation 
rates T2 and 1/t„ cause transitions in both directions with a rate (1/2t„ + T2/2)N^^i 3> FiiV, with TV = iV| + Ni. 
We go to the continuous limit, justified by the large number of nuclei per dot {N ~ 10^) to obtain [7] 



aP(i?^,S«,t)_ d f ^dB^ ^ 2Bl d I 



dt dBj^ \ dt N dBj^' V2r„ 



-r- 



dB§ , 2Bl d ^( \ 



dB§ \ dt N dB§ V2t„ 



(11) 



From the steady state solution of ( |TT| ) we evaluate the small fluctuations of the nuclear flelds around the stable states. 
For any unpolarized dot we find ((Ai?jv)'^) = = fl^, i.e. the fiuctuations are not affected by ESR. If one or both 

of the dots are polarized, then we can express the resulting nuclear field fiuctuations in the polarized dot in terms of the 
maximally reachable field as ((ABtv)^) ~ {Bi/\B'^^^\)^l'^, i.e. the fiuctuations are suppressed by a factor Bi/lB'^'^^]. 



MORE DETAILED RESULTS 



Here we will present three plots in addition to Fig. 4 in the main text. The plots in this section are generated using 
the same parameters as in Fig. 4, but supply some extra details which were omitted from Fig. 4 for reasons of clarity: 
Here we include vector field plots of the time-derivatives {Bj^,B^} in the plane {Bj^,B^). Starting from a specific 
nuclear field configuration {B^, B^), following the arrows shows the evolution in time of the nuclear fields. We added 
the current through the system /dot as color background, this gives a more quantitative picture of the current levels 
in the different stable points. 

As mentioned above, the current through the double dot can be calculated from the quasi-stationary populations 
Pn as /dot = eFs ^„ |(S'|n)pp„, and is a function of fL,R/ B^ji (and therefore of the nuclear polarizations in the two 
dots), the temperature S, and the ratio Tr\Ez\lTs- In Fig.lTlwe plotted /dot close to the point where both dots are on 




FIG. 9: In the middle of the hysteresis interval, Bo ~ Bics + 0.5 assuming a symmetric double dot (Fig. 4d in the main 

text). Again we show vector plots of {B^,B^} in the plane {B^,B§). (a) Overview of the whole plane, where the circles 
indicate the stable points of nuclear polarization. Owing to the double-peak structure of the pumping curve, the left(right) dot 
has three stable points along the line B^'^' = 0, i.e. where the right(left) dot is unpolarized. (b) Close up of the region where 
the right dot is unpolarized and the left dot is close to resonance. One of the stable points in this region corresponds to a high 
current through the system, the other to low current, (c) Close up of the region where both dots are close to resonance. Four 
additional stable points can be distinguished, two of which correspond to low current and two to high current. To generate 
these plots, we used ^rr[Ez]/r, = 0.75, Fi/Fa = 21 ■ 10-^ FaTn = 5, and Bo - cu = 0.475]^''''. 



resonance, i.e. where fL = fR = 0. The function has the structure of two crossing Lorentzians, with a suppression at 
the resonant point /l = /i? = 0. In all current plots in this section, we subtracted the leakage current far away from 
resonance: It is a measure for the spin relaxation rate rj.[i?z]- 

In Fig.jsjwe replotted Fig. 4a from the main text, and added {S^, B§} as vector field and the current /dot as color 
background. In this case Bq — Bios ^ Bi, i.e. the detuning of Bq and uj is relatively small. In the whole plane we 
distinguish four stable points: two with low and two with high current. As can be seen from the close ups in Figs[8]3 
and c, the two values of high current can differ with ca. 15 %. The asymmetry in B^ f; and N^^n is implemented by 
using rf/rf — 0.097 and r2 /Ff = 0.44. This corresponds to a difference in a.c. magnetic fields B^'^ of ~ 50 % and 
a difference in effective numbers of nuclei of ^ 30 %. 

In Fig.|9]we present the same plot as in Fig. 4d in the main text, again with the vector field of time evolution and 
the current added. Here Bq « Bres + 0.5 i.e. the system is in the middle of the hysteresis interval and we 

assume a symmetric double dot, i.e. equal parameters for both dots. In the whole plane nine stable points can be 
distinguished. In four of those points the current through the double dot is relatively high. The unpolarized point 
B^ = B^ = is so far away from the other stable states that, as soon as the system switches to the unpolarized 
state, it will stay there forever. 



7 



[1] F. H. L. Koppcns, ct al., Nature 442, 766 (2006). 

[2] R. Hanson, ct al. Rev. of Mod. Phys. 79, 1217 (2007). 

[3] J. Danon and Yu. V. Nazarov, Phys. Rev. Lett. 100, 056603 (2008). 

[4] A. W. Overhauser, Phys. Rev. 92, 411 (1953) 

[5] M. S. Rudner and L. S. Levitov, Phys. Rev. Lett. 99, 246602 (2007). 

[6] In this expression we assume for simplicity nuclei with spin 1/2. In GaAs I = 3/2, which gives rise to an extra numerical 

prcfactor in the later results. 

[7] N. G. van Kampen, Stochastic Processes in Physics and Chemistry (North-Holland, Amsterdam, 1990). 



